Macroscopic thermal entanglement due to radiation pressure 
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Can entanglement and the quantum behavior in physical systems survive at arbitrary high tem- 
peratures? In this Letter we show that this is the case for a electromagnetic field mode in an optical 
cavity with a movable mirror in a thermal state. We also identify two different dynamical regimes 
of generation of entanglement separated by a critical coupling strength. 

PACS numbers: 



Recently there has been a lot of interest in macroscopic 
systems which can support non zero temperature entan- 
glement (l| |2j LI • These works suggest that entanglement 
not only persists in these conditions but that it is also 
crucial in understanding the macroscopic proprieties of 
the physical systems. On the other hand there is also 
a possibility that, under some circumstances (very large 
number of atoms, very large temperature, etc.), the phys- 
ical systems become effectively classical and we no longer 
require a quantum description. 

It is usually believed that decoherence explains the 
transition from the quantum world to the classical world. 
Several other alternatives, such as spontaneous wave- 
function collapse models |j| and gravitationally induced 
decoherence ja] are also found in the literature. Differ- 
ent schemes to create and probe macroscopical superposi- 
tions, which can give us important clues about the quan- 
tum to classical transition, have already been proposed 
for different physical systems 00 • For instance, in 
, like in the original gedanken experiment proposed by 
Schrodinger, a single photon state induces quantum su- 
perpositions of a mirror and in pj multi-component cats 
of a cavity are created in the interaction of a cavity field 
with a movable mirror. 

Here we proceed in a different direction shedding light 
into the problem of how quantum features may survive 
in the macroscopic world. Recent research on this prob- 
lem have shown that in systems with finite Hilbert space 
dimension in equilibrium with a thermal bath, bipartite 
entanglement vanishes above a critical temperature |lfj| . 
The following question arises, is this behavior general, 
or are there some systems where entanglement is robust 
against temperature? 

In this Letter we will show that the entanglement be- 
tween macroscopic mirror and a cavity mode field can 
arise due to radiation pressure at arbitrarily high tem- 
peratures as the system evolves in time. This is very sur- 
prising because it is commonly believed that high temper- 
ature completely destroys entanglement. We will study 
entanglement in the time domain using a discrete vari- 
able method and identify its dependence on the relevant 
physical parameters, such as the strength of the radiation 



pressure coupling and temperature. 

In particular, we consider a perfect optical cavity with 
a movable mirror with mass m in one end, modeled as a 
mechanical harmonic oscillator whose quivering motion is 
quantized. Identical systems have been considered in the 
studying of decoherence [5| and non classical states of the 
cavity field . The general Hamiltonian of this system 
has been extensively studied by Law 0] . Following the 
same approach as Law, we consider the adiabatic limit 
where the resonant frequency of the mirror is much slower 
than the frequency of the cavity mode w m << 2imc/L, 
where L is the length of the cavity when the mirror is 
in equilibrium, n is the order of the longitudinal cavity 
mode and c is the speed of light. Henceforth, the cou- 
pling between different cavity field modes (leading to the 
Casimir effect, etc.) can be neglected. For cavities with 
very high quality factor Q the damping is negligible, as 
it occurs on a time scale much longer than it takes for 
the photons to perform several round trips. Under these 
conditions the Hamiltonian includes only the free terms 
of both the field and the mirror plus the interaction term 
due to the radiation pressure (which causes the displace- 
ment of the mirror), 



H = hw$a) a + hw m wb — hga) a(b + b^), 



(1) 



where a is the annihilation operator of the cavity mode, 
b is the annihilation operator of the mirror, wq is the 
frequency of the cavity mode and g — woVh/ (Ly/mw m ) 
is the coupling constant. 

The electromagnetic field is prepared in a coherent 
state, \a), using a driving laser tuned to resonance with 
the cavity mode, whereas the mirror is initially in a Gibbs 
state with temperature T. Then, the composite state of 
the system is 
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\a)(a\ ® \z)(z\, 



where n = \j{e hw ^l K BT _ ^ j g mean number of 
excitations, Z is the mirror partition function and z rep- 
resents all the possible coherent states of the mirror. 

The evolution operator associated with the Hamilto- 
nian (1) has a closed formula and it was derived in 0, 
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using the Campbell-Baker-Hausdorff formula for the Lie 
algebra, and in J\, using operator algebra methods: 



U(t) = e - lWoatat e' k2 ^ a '> 2A ^D m (r 1 (t)k^a)e- lWmb1bt , 



where A(t) = w m t — sin(w m t), r)(t) 
k = g/w m and D m (r)(t)ka^a) = e ^ a( n (tp -r,(t)'b) ig the 
displacement operator of the mirror, D m {-~{)\0) = \ r y). 
Since the system is periodic we only need to investigate 
entanglement in the time interval [0, 27r/w; rra ]. 

The interaction term of the Hamiltonian has the po- 
tential to entangle the cavity field modes with the vibra- 
tional modes of the mirror. Heuristically, the generation 
of entanglement in this physical system is better under- 
stood by considering the cavity in the initial state |0) + 1 1) 
and the mirror in the vacuum state |0) m (which is a good 
assumption for T « 0). The state of the composite sys- 
tem evolves according to (up to a normalization factor) 

(|0) + 11» <8> |0) m -> |0) <g> |0) m + e^|l) ® \kr,(t)) m , 



where f(t) is a phase, resulting in an entangled state for 
< t < 2n/w m . At t = and = 2it/w m , 7](t) becomes 
null and the system returns to the initial separable state. 
The entanglement results from the evolution of the term 
|1)<8> |0) m , which can be interpreted as the transference of 
momentum from the photon |1) to the mirror |0) m , as the 
photon kicks the mirror. Though it is easy to produce the 
state of light |0) + |1) experimentally, the radiation pres- 
sure in this case is so small that its virtually impossible 
to detect any entanglement using present day technology. 
However, as we will show below, a detectable amount of 
entanglement is expected when the cavity is initially in 
a coherent state with sufficiently high amplitude. 

Proceeding to the general case, we introduce the dis- 
crete variable method which will allow the study of the 
entanglement in the system. The density matrix p(t) ex- 
pressed in the Fock basis of both cavity field and mirror 
reads: 



^ Pnunm\n)(m\ ® \fJ,)(v\ 



(2) 



where the Latin indexes refer to the radiation and the 
Greek indexes refer to the mirror. The elements of the 
density matrix are calculated to be: 



, for t = 2 / nk/w m or n = m = 



P^ 



W^Jty-KHui-n, 1 + z nm (t)} 
, elsewhere. 



(3) 



Hjj[a 1 b, c] is the hypergeometric confluent function and 

a n a *m e iA(t)(n 2 -m 2 )-\a\ 2 

W nm (t) = A n (t) - e-^iAnit) + A m (t))/2 
z nm (t) = -e 0hw ™W nm {t)W^ n {t) 
n nm (t) = (|A„(t)| 2 + |A m (i)| 2 )/2 

_e-^-|A n (t) + A m (t)| 2 /4 

where A n (t) = knr)(t). 

Quantifying entanglement in mixed states is a non triv- 
ial problem, except for bipartite two level systems, where 
the Peres criterium is sufficient and necessary to 
guarantee entanglement. In this Letter a new approach 
inspired by 0| is followed. First, we project the original 
density matrix (2) into a 2 x 2 subspace, which corre- 
sponds to a local action (P2X2P-P2X2) thus not increasing 
the amount of entanglement E(p) in the overall system 
[3, i.e. E(p) > E(P 2x2 pP2x2)- Afterwards different 
measures and markers of entanglement for 2x2 systems 
can be applied. In particular, the tangle is an entan- 
glement monotone valid only for bipartite pure states 
|l7|. The negativity is also an entanglement monotone 
but valid for both pure and mixed bipartite states 0] 
N(p) — Q2i M ~ l)/2, where are the normalized 
eigenvalues of the partial transposed matrix p Tp . The 
relation between negativity and the Peres criterium is 
clear: it measures by how much the partial transposed 
density matrix fails to be positive, i.e. separable. There- 
fore to show the existence of entanglement in the overall 
system it is sufficient to verify that E(P2 X 2pP2x2) > 0. 
From now on we will use the notation where each sub- 
space of the density matrix is represented by [p, v\ n, m] 
where p, v refer to the number of excitations of the mir- 
ror and n, m refer to the number of excitations of the 
cavity field. 

For T = the system is in a pure state and we can in- 
vestigate entanglement using the tangle [13, r(t), which 
is much simpler to compute than the negativity. Fig- 
ure 1 shows that the system is always entangled except 
for t = and t = 2ixjw m (only the range t € [0, ir/w m ] 
has been plotted since the tangle has reflection symmetry 
around t = ■n/w m ), when the mirror returns to its ini- 
tial state. For small k, the system reaches the maximum 
of entanglement at t = ■n/w mi simultaneously with the 
maximum displacement of the mirror. For k above a crit- 
ical value, k c , the maximum of entanglement is achieved 
before t = n/w m . The time of maximum entanglement 
depends on the balance between the interaction time i,- n < 
(tint ~ 1/5)1 i- e - the time scale of the interaction term 
in the Hamiltonian, and the time of oscillation of the 
mirror, t m (t m ~ 1/wm). 

It is interesting to try to understand the importance 
of the amplitude of the coherent state, a, in establish- 
ing the value of k c . We expect that when increasing a 
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the value of k c should decrease because there are more 
photons interacting with the mirror resulting in a larger 
effective coupling (g((va)). Surprisingly this is not the 
case: the value of k c increases with a. This can be un- 
derstood as follows. The ratio between the weight of the 
\n + 1) number state and the weight of \n) number state 
in the expansion of the coherent state, being given by 
a/{s/n + 1), increases with a, weakening the entangle- 
ment generated after interaction with the mirror. The 
best situation occurs when the weights of states are the 
most equally distributed. Hence a higher coupling helps 
the entanglement generation to have the same efficiency 
when a is increased. 

Considering the subspace [1,2; 1,2], k c can be calcu- 
lated as solution of the transcendental equation 

- 1 + Uk 2 c + 24kj = 2a 2 (l + 4fc 2 + 96kf) e - 12k * > 0.(4) 

The right-hand side of equation (4) is non-negative re- 
sulting in a restriction for fc c , i.e. k c is lower bounded. 
Also we can see from equation (4) that k c increases with 
a. This confirms, at least for this subspace, that a higher 
coupling is necessary for reaching the maximum of entan- 
glement before t — ir/w m if the amplitude of the cavity 
field is increased. 

Since t(tt) is the maximum of entanglement achieved 
during the evolution of the system for k < k c , and k c 
depends on a, the value of t(tt) is maximized for 

_ e 6fc Vl + 2fc 2 
amax ~ V2Vl + 8fc2 ' (5) 

which can only be defined for k c {a max ) > k. Equation 
(4) indicates that the ideal value of a increases with the 
coupling. It is useful to define a c as the value of a, for a 
given k, such that k c — k and to notice that k c {a max ) > k 
is equivalent to the condition a c > a max . By squaring 
(4) and dividing it by (5) (with k c — > k and a — > a c ) the 
condition of validity of (5) yields k > 1/2. This behavior 
is independent of the particular subspace considered. For 
example, the asymptotic behavior of the tangle at t = n 
as function of a is always t(tt) ~ |a| 2 , for a <§; a max , 
and t(tt) ~ | | 2 , for a 3> ot max . Figure 1 shows how 
the maximum of entanglement in the considered subspace 
is shifted to earlier times when k increases. 

For temperatures slightly above T = this behavior 
is not significantly altered. At temperatures T > the 
system is in a mixed state and the entanglement must be 
investigated using the negativity. It can be inferred from 
the plots of the negativity obtained from the simulations 
that, for the subspace [1, 2; 1, 2], the value of k c increases 
slowly with the temperature. The investigation of the 
negativity also shows, as expected, that the entanglement 
decreases with temperature in a given subspace. 

For higher temperatures, i.e. higher number of excita- 
tions of the mirror, the negativity cannot be computed 




FIG. 1: Tangle r as a function of k and the scaled time for 
T = OK, a = 1 and subspace [1, 2; 1, 2]. 

because the numerical calculation of the hypergeomet- 
ric confluent function is ill conditioned. However it is 
still possible to verify the existence of entanglement for 
different values of a, k and temperature. This can be 
accomplished by introducing a marker of entanglement 
based on the Peres criterium, which consists simply in 
verifying whether the determinant of the partial trans- 
posed density matrix p Tp is negative, thus indicating the 
existence of a negative eigenvalue. This is of course a 
very weak marker of entanglement, which can only de- 
tect entanglement if p Tp has a odd number of negative 
eigenvalues, but has the useful advantage of being easy 
to calculate. 

The determinant of p Tp for the subspace [0, 1; 0, 1] is 
T[o ) i ; o,i](q:, b, x) = —G(a 1 b,x)H(b,x), where H(b,x) — 

16x 2 e b2x + x 2 (-4 + b 2 (-2 + x){-l + e 6 '/ 2 )) 2 - 
e 6 2 */ 2 ( 6 4(_ 2 + x) A + 32.x 2 + 46 2 a; 2 (-l + e fc2 / 2 )(4 + (-3 + 
e b2 / 2 )x)), G{a,b lX ) = x 4 |q| 4 e b2 ( x - 2 )- 4 H 2 /(16 7T 4 ) is al- 
ways positive, b = b(t) — \[2k\J\ — cos(w m t) and x = 
71/ (1 + n). The sign of T is determined by H (b, x) and is 
plotted as a function of (2 / ' -n)arctan(b(t)) and x in Figure 
2, where it is clear that the entanglement occurs for all 
parameters except for high values of b(t) and low tem- 
perature (region I) and for low values of b(t) and high 
temperatures (region II). 

In the parameter region I, the lack of entanglement is 
not important since, during the evolution of the system, 
the value of b(t) ranges between and b max — 2k and, for 
sufficiently short times, the system is entangled indepen- 
dently of how large k is. In the parameter region II the 
proof of existence of entanglement in the system is not 
as easy because, in principle, the system could always be 
separable for small values of k. However, by rewriting 
equation (3) as follows 

J I" 

where K nm {z) = -(|F„(z)| 2 + \F m (z)\ 2 )/2 - \z\ 2 /n and 
F n (z) = z + knr](t), it is clear that k is always multiplied 
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either by n or by m. Then it is straightforward to verify 
that TCVijQ gi it is proportional to T[ 0i i ; o,i] if in the latter 
the coupling k a = k/s is chosen, i.e. 

|ar 2s S !T [0!l;0 , s] (fc s ) = |a|- 2 T [0!l;04] (fc). 

For high temperatures, if we choose a large k capable of 
producing entanglement in the subspace [0,1; 0,1], then 
there must be entanglement in the subspace [0, 1; 0, s] for 
the coupling constant k s = k/s, even though k s might not 
lead to entanglement in the subspace [0,1; 0,1]. There- 
fore, for systems where b(t) varies in region II there is 
always a s such that entanglement occurs in the subspace 
[0,1; 0, 5 ]. 




0.2 0.4 0.6 0.8 1 
(2 / it) arctan (b) 



FIG. 2: Marker of entanglement, T, as a function of 
(2/tt) arctan(6(t)) and x for the subspace [0,1; 0,1]. Entan- 
glement (gray) exists for all parameter space except in regions 
I and II (white). 

Although entanglement occurs at any finite tem- 
perature this does not imply the existence of en- 
tanglement in the limit of infinite temperature. 
In fact, making use of the general inequality 
E(Y,iPiP\ B ) < EiPiE(p\ B ) to write E{Pp{t)P) < 
J dze-W 2 / n E{PU\a)(a\®\z)(z\WP)/nTT and taking 
into account that i) the limit of infinite temperature is 
equivalent to the limit n — > oo and ii) entanglement is 
bounded in any subspace, yields Yimn^oo E(Pp(t)P) = 0. 

The entanglement appears naturally in the dynamics of 
the system for any value of coupling between the mirror 
and the cavity field at any finite temperature, indepen- 
dently of the value of a (notice that H does not depend 
on a, at least for the class of subspaces considered here), 
in spite of a playing an important role in determining 
the amount of entanglement in the system. An argu- 
ment will be given to infer this role. We define the nor- 
malized mutual information as X = lyxj [S\ + S2), where 
/12 = Si + S*2 — S, Si is the Linear Entropy of the mirror 
(i = 1) and the cavity field (i = 2). 5 = 1- l/(2n+ 1) 
is the Linear Entropy of the composite system and 

oo I \2(r)-\-q) 

St = l- e- 2 ^Mn) Y H^ e - Sl (a, fc ,w, t )(p- 9 ) 2 ; 



for t ^ 0,27r/w m , where /i(n) = l/(n 2 (a 2 - 1)), / 2 (n) = 
1, gi{a,k,n,t) = 2fc 2 £(i)(l - b), b = (1 - l/af + l/a, 
a = 1 + 1/rT, = 1 — cos(iu TO t) and g 2 (a, k,n,t) — 
k 2 £(t)(n + 2). Since gi(a, k,n,t) > (for k ^ 0) we have 
dl/d\a\ > 0, which shows that quantum correlations in- 
crease monotonously with a as mentioned before. Hence, 
a detectable amount of entanglement is expected for suf- 
ficiently high a. Finally, we give a numerical value for I 
for realistic physical parameters T = lK,w m = 10MHz, 
k = 1 and a = 10 6 ,X-l/2 > 0.19xl0~ 4 > 0. TheAraki- 
Lieb inequality 18] guarantees that quantum correlations 
are present whenever I > 1/2. This is consistent with 
our previous conclusions that entanglement persists even 
in the high temperature limit. 

In this Letter the generation of entanglement between 
a cavity mode and a movable mirror via radiation pres- 
sure was analyzed by a projection method. Since in the 
process of projection of the density matrix much of the 
entanglement is lost, this method does not allow to assess 
the entanglement of the overall system but it is still ro- 
bust enough to identify two dynamical regimes depending 
on the coupling between the mirror and the cavity field. 
For k < k c the maximum of entanglement is achieved 
at t — 7r and for k > k c is achieved before and further- 
more demonstrate the remarkable result that entangle- 
ment does occur for any finite temperature. 

We hope that our results stimulate further investiga- 
tions on entanglement in the macroscopic domain. We 
gratefully acknowledge M. Aspelmeyer, C. Brukner, S. 
Gigan, J. Kofler, E. Lage, P. Vieira and M. Wiesniak. 
A.F. is supported by FCT (Portugal) through grant 
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